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ABSTRACT 

A  Hugoniot  pressure  density  relationship  for  homogeneous  mixtures 
la  calculated  from  the  Hugoniot  pressure  density  relationship  of  the 
constituents.  The  c< Iculated  values  are  found  to  be  in  reasonable 
agreement  with  available  experimental  data  on  mixtures  of  solids  and 
porous  materials. 
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I.  INTRODUCTION 

The  motivation  for  the  study  here  reported  is  the  intersection 
of  two  trends  in  technplogy;  one  being  the  increasingly  wide  appli¬ 
cations  of  new  composite  materials;  the  other  being  the  rapid  growth 
of  interest  in  the  ^response  of  structures  and  structural  elements  to 
extremely  strong  disturbances. 

Shock  properties  of  mixtures  have  been  computed  [1]  from  the 
properties  of  the  constituents  by  a  method  using  the  complete  equation 
of  state.  In  this  report,  a  means  of  computing  approximate  shock  pro¬ 
perties  from  the  Uugonlot  relationships  alone  is  proposed.  The  results 
of  this  study  will  be  compered  with  experimental  data  for  an  unstruc¬ 
tured  or  homogeneous  composite.  The  simplified  theory  will  also  be 
investigated  as  a  means  of  predicting  the  shock  response  of  porous 
materials. 

In  an  earlier  study,  a  theory  for  predicting  the  response  of  a 
layered  composite  to  a  propagating  shock  was  given  [2].  It  is  to  be 
expected  that  the  particular  geometrical  arrangement  of  the  two  consti¬ 
tuents  assumed  in  that  study  would  strongly  Influence  the  propagation. 

In  the  present  study,  the  propagation  of  a  shock  wave  through  a  compo¬ 
site  material  having  no  geometrical  structure  is  of  Interest.  It  is 
anticipated  that  a  comparison  of  the  results  from  these  two  theories 
will  provide  at  least  a  qualitative  Indication  of  the  influence  of 
structure  on  the  propagation. 

A  related  simple  theory  for  the  speed  of  propagation  of  an 
acoustic  disturbance  is  also  considered  and  Is  found  to  display  certain 
interesting  features.  The  condition  under  which  the  speed  of  sound  in  a 
mixture  can  be  leas  than  the  sonic  velocity  of  either  constituent  is 
established. 


II.  THE  HUGONIG!  RELAX iUNSHiP 

We  wish  Co  determine  the  shock  response  of  a  homogeneous  composite, 
such  as  would  result  if  powders  of  several  constituents  are  prepared, 
thoroughly  mixed,  and  then  pressed  so  as  to  fora  a  solid  having  the 
desired  mass  (or  volume)  fractions  of  the  constituents.  The  resulting 
composite  is  assumed  to  be  macroBCopically  homogeneous,  l.e.,  it  behsvss 
as  a  homogeneous  and  isotropic  solid  and,  to  the  scale  of  the  objects 
which  might  be  made  from  the  composite  material,  no  structure  can  be 
seen.  This  criterion,  it  should  be  noted,  may  not  be  met  In  a  shock 
process,  for  the  distance  over  which  the  load  ia  applied  is  the  thickness 
of  the  shock,  the  product  of  the  rise  time  of  the  shock  and  the  speed 
of  propagation. 

a.  Theory  for  a  Mixture  of  Two  Constituents 

We  assume  that  any  element  of  volume  AVo  at  one  atmosphere  of 
pressure  contains  a  volume  AVj,  of  one  constituent  and  AV2  of  the 
other.  If  the  two  constituents  have  initial  densities  p«  and  po2  , 
then  the  masses  of  each  of  the  two  constituents  are 


Mj  ■«  AVj  and 

M2  -  AV 2  poj 

The  initial  average  density  of  the  composite  is 


(1) 


ave 


Mi  ±  M?  .  ,  , 

AV  +  V#a 


(2) 


or 


(3) 
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•eve 

where  fr2  and  fVj  are  initial  volume  fractions  of  the  two  components, 
and  fmi  and  fm2  are  the  mass  fractions.  It  is  implicitly  assumed 
that  the  mixture  is  simple,  or  non-reacting. 


k 


The  average  density  of  the  composite  after  the  passing  of  the 
ahnrU  ia 


M  i  +  M  2 

Pfava  =  W\i  +  &  £ 


(A) 


where  M  t 
and  LV  jf 
Since  Mj 
the  shock. 


and  M2  are  the  masses  of  each  constituent  and  AV  ^ 
the  volumes  occupied  by  these  masses  at  the  shock  pressure, 
and  the  mass  fractions  fm^  are  assumed  to  be  unchanged  by 
the  possibility  of  phase  change  under  pressure  is  excluded. 


Thus 


_ _  fml  . 

Pfave  W|  ^2 


(5) 


Let  us  now  assume  that  the  response  of  each  constituent  to  a  one 
dimensional  shock  (the  Hugoniot  or  shock  adiabat)  is  known.  For  moder¬ 
ate  pressures  (hundreds  of  kilobars)  this  relationship  can  be  approxi¬ 
mated  for  many  materials  by  the  expression  [3] 

P  -  Aj(p/p0^n-lj  (6) 

where  P  is  the  pressure  behind  the  shock,  measured  above  one  atmosphere, 

A  and  n  are  parameters  of  the  material 
p  is  the  density  behind  the  shock  and 
p0  is  the  density  in  front  of  the  shock. 

A  relationship  of  this  form  is  uoually  satisfactory  only  if  no  phase 
changes  occur  over  the  pressure  range  where  it  is  to  be  applied.  Assuming 
that  the  density  of  each  constituent  of  the  mixture  at  a  given  shock 
pressure  is  the  same  ea  would  be  found  in  a  homogeneous  sample  at  the 
same  pressure,  end  that  each  constituent  can  be  described  by  Equation  (6), 


where  the  pressure,  P,  is  the  same  in  each  constituent,  and  the 
are  densities  behind  the  shock.  Substituting  (7a)  and  (7b) 
and  (3)  into  (5) 


which  provides  a  means  of  predicting  the  average  density  vs  pressure 
relationship  for  a  homogeneous  mixture  of  two  constituents,  and 
requires  only  the  knowledge  of  an  approximate  Hugoniot  relationship 
for  each. 

The  assumed  pressure  density  relationship ,  Equation  (6),  was 
selected  only  as  a  matter  of  convenience.  Any  form  giving  Hugoniot 
pressure  as  a  function  of  density  which  can  be  solved  explicitly  for 
density  as  a  function  of  pressure  could  as  well  be  used  and  would  lead 
to  a  closed  form  expression  analogous  to  Equation  (8). 

For  many  materials,  Hugoniot  data  is  readily  available  (4)  In  the 

form 

D  ■  C  +  SU  (9) 

where 

D  is  the  shock  speed 
U  the  particle  velocity,  and 
C  and  S  are  empirical  constants 
It  has  been  demonstrated  that  such  data  [5]  may  be  put  in  the  form 
of  Equation  (4)  through  setting 

Poc2 

A  *  ;  -  and  n  -  4S-1 

4S-1 


(10) 


It  iias  alsv  been  shown  f 3 j  cnac  universal"  values  of  A  and 
n  may  be  used  with  good  accuracy  for  a  wide  variety  of  materials. 

These  values  are  n  ■  5  and  A  *  p0C02/5.5,  where  is  the  initial 

density  and  C#  is  the  sonic  velocity  at  standard  conditions.  By 
using  these  properties  of  the  constituents.  Equation  (8)  can  be  used  to 
predict  the  Hugoniot  pressure  volume  relationship  even  if  Kugonloc 
information  for  the  constituents  is  not  available. 


b.  Response  of  a  Hypothetical  Copper-Polyethelene  Mixture 

The  family  of  predicted  shock  pressure  vs  density  (  '>  relation¬ 
ships  for  a  hypothetical  copper-polyethelene  composite  was  determined 
through  Equation  (8).  A  linear  Hugoniot  relationship 

D  -  3.92  +  1.488  J  <U) 

with  P0  -  8.93  gm/cm  was  chosen  [6]  for  copper.  All  velocities  are 
in  units  of  km/sec._  Equation  (10)  then  yields  parameters  A  •  279.7  kb 
and  n  ■  4,956.  Available  data  [1]  on  polyethelene  was  used  to  deter¬ 
mine  the  coefficients  to  be  A  ■  9. 64  kb  and  n  ■  5.875.  The  initial 
density  was  taken  to  be  .915  gm/cm  .  The  results  of  these  calculations 
are  shown  in  Figure  1  as  pressure  vs  density  relationships  for  various 
mass  fractions  of  copper.  A  mass  fraction  of  .9,  for  these  materials, 
corresponds  to  a  mixture  of  48Z  copper,  by  volume. 

Once  the  Hugoniot  pressure  vs.  density  relationship  has  been  computed, 
the  other  shock  properties  can  be  determined  from  the  Rankine-Hugoniot 
jump  conditions,  i.e. 

Dp0  -  (D-U)Pf  (12a) 

r°ave  rave 

P  »  Po  DU  (12b) 

ave 

where  the  composite  is  assumed  to  be  macroscopically  homogeneous.  Qi.e 
particularly  interesting  result  is  given  in  Figure  2,  where  the  shock 
speed  at  various  pressure  levels  is  plotted  against  the  mass  fraction  of 
copper.  The  pronounced  minimum  at  low  pressures  is  particularly  signifi¬ 
cant  and  suggests  the  possibility  of  a  minimum  in  the  sonic  velocity. 

This  subject  will  be  treated  in  a  later  section. 

c.  Comparison  with  Theory  for  a  Layered  Composite 

In  Figure  3,  shock  speed  vs.  particle  velocities  as  computed  from 
this  theory  (the  solid  lines)  are  compared  with  results  obtained  previously 
for  a  layered  composite  having  the  same  composition.  Several  mixtures 


of  aluminum  with  an  assumed  linear  Hugoniot 


<;  o/. 


and  polymethylmethacrylate ,  with  an  assumed  linear  Hugoniot  of 
D  -  2.70  +  1.61U  km/sec 


were  considered.  The  results  showed  a  surprisingly  good  agreement 
between  the  two  computations,  suggesting  that  the  influence  of  the 
layering  assumed  in  [2]  is  much  less  than  might  be  expected. 

d.  Theory  for  Several  Constituents 

The  theory  may  be  readily  extended  to  mixtures  of  three  or  more 
constituents.  Denoting  the  1  th  constituent  by  a  subscript  i,  the 
mass  of  that  constituent  in  some  volume  element  AV  is 


Then 


and 


Mi  “  Avi  P«i  "  WVP  1 


“avo 


ar-LIll. 
E  AV, 


.  ? 


P>L 


m _ 1 
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(13) 

(14) 


1-1  P®i 


^“i  '  ^V-I 

Poave 


Assuming  each  constituent  may  be  described  by  Equation  (6)  or 


since 


_1  _  ^tm 

.  f£i 


ave 


(7c) 


III.  APPLICATION  OP  THEORY  TO  POROUS  MEDIA 
Of  particular  Interest  are  composites  where  one  of  the  constituents 
la  air.  such  as  would  result  from  an  Imperfect  packing.  For  air,  hov= 
ever.  Equation  (6), is  not  appropriate.  We  assume  air  to  he  a  perfect 
gas  with  constant  ratio  of  specific  heats  y,  end 

P  -  pRT  (15a) 

E-$>  -  Cv(T-T^)  (15b) 

Using  these  and  the  Ranklne  Hugonlot  Jump  condition  for  a  one  dimen¬ 
sional  shock  in  a  homogeneous  material, 

E-Eb  ■  |Vp0  -  Vpfj-  (l6) 

to  determine  a  shock  pressure  vs  density  relationship  for  air  which  is 
analagous  to  Equation  (6),  we  find,  for  y  *  1»^ 

|  -  ^Oo  '  1  or  g.  -  (17) 

P°  6  -  P/p0  Po  6  +  P/po 

Thus  the  average  density  of  a  composite  containing  N  solid  constituents 
and  air  is 


N  D  -1 
K®ave  v-  -  rP  .  i  - 

i  vifc  /  1  v« 

ttV*  i-1 


r  6  +  p/p0 

&  P/P0  + 


where  fva  is  the  initial  volume  fraction  of  air.  For  pressures  of 
interest  for  solids,  P/P0  »  1  and  the  last  term  may  be  written  as 

*^a/6.  Experimental  data  [7]  on  strong  shocks  in  air  indicates  that 
the  compression  under  shocks  of  strength  greater  them  200  bars  is 
closer  to  ten  than  the  value  of  six  predicted  by  treating  air  as  a 
perfect  gas. 


a,  Comparison  of  Theory  with  Experimental  Results  on  Materials  of 
Low  Porosity 

Experimental  data  for  a  tungsten-copper  mixture  (EUccalte)  was 
obtained  from  the  literature  [i]  and  compared  with  values  predicted 
by  this  theory  using 

D  -  (4.029  +  1.237U)  km/sec  (19) 

with  pG  •  19.22  gm/cm3  for  tungsten  and 

D  -  (3.940  +  1.489U)  km/sec  (20) 

with  p0  -  0.93  gm/citf  for  copper, 

Constants  A  »  .790  mb,  and  n  -  3.948  were  determined  for 
tungsten  and  constants  A  ■  .280  mb  and  n  -  4.956  were  determined  for 
copper. 

The  average  compositions  of  the  mixtures  for  which  Qugonlot  data 
was  available  are  given  in  Table  1.  Shock  pressures  corresponding  tr> 
the  theoretical  densities  given  in  Table  I  were  computed  and  compared 
with  the  experimental  data.  The  agreement  between  the  experimental 
data  and  the  elementary  theory  was  found  to  be  remarkably  good.  It 
can  be  seen,  however,  from  Table  I  that  the  samples  apparently  have 
some  porosity;  that  Is,  the  measured  densities  are  not  precisely 
the  theoretical  density  of  a  mixture  having  the  stated  mass  fraction 
of  the  two  constituents,  but  are  always  somewhat  lower.  If  the  differ¬ 
ences  can  be  attributed  to  porosity,  and  the  expression  given  by 
Equation  (18)  takes  proper  account  of  the  pressure  of  air,  then  Equation 
(10)  should  yield  values  of  final  density  which  are  in  better  agreement 
with  the  experimentally  determined  values  than  predictions  as  based  on 
Equation  (8) .  For  each  composition  (nominal  mass  fraction  of  each)  the 
quantity 


prediction  of  the  final  density  may  be  obtained  by  neglecting  the 
porosity  and  using  the  initial  density  the  density  vhlch  the  com- 
puttile  would  have  if  there  were  no  porosity.  The  three  methods 
of  calculation  are  compared  in  Figure  U.  The  solid  lines  Indicate 
typical  pressure-density  relationships  as  might  be  computed  from 
Equation  (8)  if  the  measured  density  p0  or  the  theoretical  density 
(for  no  porosity)  were  used.  The  dashed  curve  shows  the  result  to 
be  expected  If  Equation  (l8) ,  taking  account  of  porosity,  were  to 
be  used.  Above  the  pressure  at  which  the  air  is  essentially  com¬ 
pressed  to  its  limiting  density,  the  difference  between  the  two 
curves  is  very  small. 

The  comparison  of  theory  using  Equation  (l8)  and  experiment 
is  shown  for  each  data  point  below  2mb  in  Figure  5*  The  predicted 
final  density  was  computed  from  Equation  (l8),  with  the  volume 
fraction  of  air  being  determined  for  each  specimen  from  the  stated 
nominal  composition  and  the  measured  density  of  that  specimen.  Aside 
from  the  systematic  deviation  in  the  U5/ 55  samples,  the  agreement  is 
remarkably  good,  considering  the  drastic  assumptions  included  in  the 
theory. 

As  was  indicated  earlier,  once  the  Hugonlot-  pressure  density 
relationship  is  established,  the  other  shock  parameters  may  be 
determined.  The  copper-tungsten  was  found  to  have  a  minimum  in  the 
Isobars  of  the  shock  speed  vs  mass  fraction  relatlchshlp  as  was  found 
in  the  copper-polyctheiena  mixture. 


b.  Comparison  of  Theory  with  Results  on  Materials  of  High 
Porosity 

Much  data  has  been  obtained  from  shock  experiments  on  porous 
materials.  Such  data  are  of  Interest,  for  the  use  of  porous  samples 
permits  points  on  the  equation-of -state  surface  for  the  solid  to  be 
reached  In  a  shock  experiment  /which  are  not  attainable  In  a  solid 
specimen.  The  porous  sample  may  also  be  viewed  as  a  rather  extreme 
composite  material  of  the  type  being  considered  herein*  Shock  data 
on  porous  samples  of  copper,  aluminum,  beryllium,  magnesium,  and 
uranium  were  obtained  from  the  literature  and  compared  with  the 
predictions  of  Equation  (l8).  Shock  pressures  are  plotted  against 
particle  velocity,  where  the  particle  velocity  for  the  theoretical 
curve  Is  determined  through  the  use  of  Equations  (12a)  and  (12b) . 

The  results  of  this  comparison  are  shown  in  Figures  6,  7,  and  8.  In 
all  of  these  comparisons,  fv  Is  the  volume  fraction  of  the  solid. 

The  comparison  of  theory  and  experiment  for  copper,  using  data  from 
McQueen,  et  al  [l]  shows  a  rather  good  agreement,  particularly  below 
250  kb.  Experimental  data  for  pure  copper  (fv  -  l.o)  have  been 
Included  for  purposes  of  comparison.  The  theory  generally  under¬ 
estimates  the  shock  pressure  in  a  porous  material. 

In  Figure  7>  data  from  McQueen,  et  al,  has  been  combined  with 
the  data  from  the  Shock  Wave  Compendium  [  6] .  The  agreement  here  Is 
not  as  good  as  In  the  case  of  copper,  the  theory  again  leading  to  a 
systematic  underestimation  of  the  shock  pressure.  The  data  for  pure 
aluminum  is  presented,  for  comparison,  as  the  solid  line,  limited 
data  for  beryllium,  magnesium,  and  uranium,  all  from  the  Compendium. 


are  compared  with  the  theory  in  Figure  8.  The  agreement  here  Is 
again  quite  good  for  all  three  materials.  Here  the  dashed  lines 
represent  data  for  the  solid  material,  while  the  solid  line  repre¬ 
sents  the  predictions  of  the  theory,  and  the  indicated  points  are 
data  from  the  literature. 

All  of  these  calculations  were  made  using  Y  *  1.4  (limiting 
compression  of  six)  for  air.  Repeating  the  calculations  using  a 
Uniting  compression  of  ten  for'  air  gave  an  insignificant  change  in 
the  results  when  plotted  as  pressure  vs  particle  velocity  or  pressure 
vs  shock  velocity.  The  final  density,  however,  is  affected. 


III.  THE  SPEED  OF  SOUND  IN  MIXTURES 


It  has  been  reported  that  the  speed  of  sound  in  certain  mixtures 
has  been  found  to  be  less  than  that  for  either  constituent.  In  parti¬ 
cular,  the  speed  of  sound  in  "bubbly"  water  has  been  found  to  be  as 
low  as  65  ft/sec.  [10] 

It  will  be  shown  in  what  follows  that  an  elementary  theory  for 
the  speed  of  sound  in  mixtures  can  be  developed  along  the  lines  of  the 
development  in  the  preceding  section  which  does  predict  that  certain 
mixtures  will  have  sonic  velocities  below  those  of  either  constituent. 
Moreover,  it  will  be  shown  to  be  unnecessary  to  assume  any  type  of 
deformation  not  present  in  the  acoustic  disturbance  of  the  constituents 
alone. 

We  will  restrict  our  attention  to  homogeneous  mixtures  of  two 
constituents;  that  is,  mixtures  in  which  the  particles  of  the  consti¬ 
tuents  are  substantially  smaller  than  any  other  physical  dimension  pre¬ 
sent  in  the  problem.  In  an  unbounded  medium,  the  only  dimension  present 
is  the  wavelength  of  the  acoustic  disturbance;  hence  the  theory  is  valid 
only  for  moderate  and  low  frequencies.  Further,  we  assume  the  relative 
amounts  of  the  two  constituents  to  be  the  same  at  all  pointB  and,  there 
being  no  structure,  the  mixture  to  be  isotropic. 

Let  two  constituents  of  density  p0^  and  adiabatic  bulk  modulus 
B-i  at  some  reference  pressure  and  temperature  be  mixed  so  as  to  create 
a  mixture  at  the  same  pressure  and  temperature  having  volume  fraction 
f^  of  each  constituent.  Since  the  mixture  is  assumed  to  be  homogeneous, 
any  elemental  volume,  ,  of  the  mixture  contains  s  volume 


AV^  -  f1AV0 


(22) 


of  each  constituent.  Since  the  mass  of  each  constituent  contained  in 
AVfl  is 


-  w7, 


(23) 


the  density  of  the  mixture  is  given  by 


.  Am,  +  Am, 


'ave 


AVC 


>o  f  +  P5  f 

j  1  2  2 


i.c.,  the  rule  of  mixtures. 

Under  a  pressure  change,  AP  “  P-P ,  the  volume  of  each 
stltuent  will  change,  but  not  necessarily  In  the  same  amount, 
nit Ion, 

dP  -  £i _  . 

^I1  P0i 


Hence  for  such  an  Infinitesimal  pressure  change, 


p  ‘  po  "  -l-_  (  %  -p0i) 

Poi 

The  volume  of  each  constituent  is  then 


AVi  "  ■“  flAVn 

P, 


bi 

and  the  average  density  at  pressure  P  is 

An*  +  Aire  Poave 


ave  AVj  +  AV2  f^ 


1  + 


P-Pj 


- FT 

+  — *- 


1  + 


P“P« 


Bi  *  ’  Ba 

Defining  an  effective  bulk  modulus  for  the  composite  through 


B. 


ave 


Po 


ave 

“pave  |  p«p9 

S- const 


Vie  find 


-ave 


£i  +  £2 

Bi  B2 


(24) 

con- 

By  def 1- 

(25) 

(26) 

(27) 
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from  which  it  is  evident  that  the  inverse  quantity,  the  compressibility 


tjaLiofiua  a  iulc  ui  lu1alui.cs. 

The  sonic  velocity  for  each  constituent  is  defined  to  be 


dF  i  Bi 
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(31) 


Assuming  the  propagation  through  the  mixture  to  be  adiabatic ,  a  aonlc 
velocity  can  be  defined 

o  dP  k 

cave  “  "  iaX6-  (32) 

Poave 

Substituting  the  expression  for  the  average  bulk  modulus  and  density, 
we  find 
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In  terms  of  mass  fractions,  fffll  and  f  rather  than  volume  fractions, 
the  above  results  become 
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replacing  the  modulus  with  the  sonic  velocities  of  the  constituents 
through  Equation  (32). 
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the  acoustic  impedance.  The  minimum  values  in  the  sonic  velocity 

in  such  cases  as  aluminum  and  polyethelene,  and  copper  and  polyethelene 
ara  balow  the  sonic  velocity  of  either  constituent.  However, 

other  combinations ,  such  as  copper  and  lsad^  *  do  not  display 
this  phenomenon. 

The  speed  of  sound  in  a  homogeneous  mixture  has  been  previously 
found  by  Wood  [11]  to  be 


l/c’  -  (op,  +  a  -«>e,l  (j£t  +  j^rrj  <38) 

where 

c  is  the  speed  of  sound  in  the  mixture 
a  is  the  volume  fraction  of  phase  2 

Pj>  pa  are  th®  densities  of  phases  1  and  2,  respectively, 
and 

ci,  ci  are  the  sonic  speeds  in  each  phaee,  defined  by 


Ci1  ■  where  the  derivative  is  evaluated  along  an 
appropriate  thermodynamic  path. 

If  a  minimum  in  the  sonic  speed  exists,  the  above  expression  for 
the  speed,  c,  sd  a  function  of  volume  fraction  a,  must  have  an  extremal 
value  in  the  range  0  <  tt  <  1,  The  sonic  speed  of  such  a  mixture  can 
also  be  written  as  a  function  of  the  mass  fraction,  x,  of  phase  2 
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where  the  average  density,  Pa,  of  the  mixture  is  given  by 


pa  “  ap2  +  (1  -a)pA  -  [x/p2  +  (1  -  x)/pi] 


The  sonic  speed, 
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in  terms  of  x,  is 
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where  Zi  are  rhe  {mfiAHflnreo. 

x  *  ' 

lane  co  Equation  (37).  Differentiating  Equation  (41)  with  respect  to 
X,  setting  to  zero,  and  solving  for  the  mass  traction  at  which  the 
extremal  value  occurs,  yields 


x  -  (B-2A)/BA  (42) 

where 

A  -  P,/P2  -  1  (43a) 

B  -  (Z JZJ*  -  1  (43b) 

But  the  mass  fraction  must  lie  in  the  range  Oixil.  Hence,  for  a 
minimum  in  the  sonic  speed  to  occur,  the  properties  of  the  two  materials 
must  satisfy  certain  inequalities.  These  are: 

If  A  and  B  are  both  negative,  then 

2A  B  ^  2A/ (l-A)  (44) 

If  A  and  B  are  both  positive,  Chen 


2AiBi  2A/(1-A)  for  AU  and  (45a) 

B  i  2A  for  A  >  1  (45b) 


If  A  and  B  are  of  unlike  sign,  no  extremal  value  can  occur.  Testing 
the  sign  of  the  second  derivative  of  the  expression  given  by  Equation  (41) 
we  find  that  the  extremal  values  which  occur  when  the  above  inequalities 
are  satisfied  are  minimal.  The  shaded  area  in  Figure  (^indicates  the 
range  of  values  of  density  ratios  and  impedance  ratios  of  the  two  mater¬ 
ials  for  which  a  minimum  in  the  sonic  speed  will  occur. 

This  analysis  can  be  applied  to  two  perfect  gasses,  each  satisfying 
the  adiabatic  pressure  volume  relationship 


(46) 


which  leads  to  sonic  velocities  given  by 


Ci 


(47) 


tor  two  such  gasses,  the  ratio  of  impedances  is 


1  TTpTP 
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(48) 


Thus,  for  two  gasses  having  the  same  y,  the  parameters  A  and  B 
defined  by  Equations  (43a)  and  (43b)  are  identical  and  the  sonic  velo¬ 
city  of  the  mixture  becomes 


Cave2  .  fm:  +  P°i  f9>2  (49a) 

Ci  Po2 


(49b) 


i.e.,  a  rule  of  mixtures.  It  should  be  noted,  however,  that  the  rule 
of  mixtures  applies  to  the  squares  of  the  sonic  velocities.  If,  however, 
the  two  gasses  do  not  have  the  same  Y,  then 


ip°i  ^ 

Y2Po2  ~  “  Po2 


(50) 


Thus,  if  the  denser  material  has  a  larger  Y,  a  minimum  in  the  sonic 
velocity  of  the  mixture  is  pooslble  and  will  occur  if  inequality  (45b) 
is  satisfied. 

One  other  special  case  is  of  interest.  If  any  two  materials  have 
the  same  impedance,  Z,  then  B  0  and 


Cave 
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a  rule  of  mixtures  for  the  velocities 


IV  SUMMARY 


A  simple  theory  for  predicting  the  pressure  density  relation¬ 
ship  for  a  homogeneous  composite  under  shock  loading  was  developed 
and  compared  with  some  available  experimental  data.  The  good  agree¬ 
ment  between  theory  and  experiment  suggests  that  tlr«  ■ imple  theory 
may  have  application  to  composites  of  this  type.  Since  this  theory 
which  assumes  the  composite  to  have  no  geometrical  structure  leads, 
in  the  one  case  considered,  to  results  surprisingly  similar  to  pre¬ 
diction  from  a  theory  for  laminated  composites,  it  may  be  that  the 
Influence  of  structure  is  r.ot  great  and  a  simple  theory  can  be  applied 
with  at  least  fair  accuracy  to  composites  which  do  display  a  geomet¬ 
rical  structure.  This  possibility  should  be  considered  further. 

The  new  theory  was  compared  with  experimental  data  for  highly 
porous  materials  by  treating  the  porous  material  as  a  homogeneous 
composite  of  metal  and  air.  The  comparison  shews  agreement  which  is 
qualitatively  good,  but  only  in  some  cases  are  the  predictions  entirely 
satisfactory. 

Conditions  under  which  the  sonic  speed  of  a  mixture  can  be  leas 
than  that  of  either  constituent  are  developed  from  the  well-known 
expression  for  the  sonic  speed.  For  a  non-reacting  mixture  of  two 
phases,  such  a  minimum  can  occur  if  certain  inequalities  relating  the 
ratio  of  densities  and  the  ratio  of  impedances  are  satisfied.  It  la 
necessary,  but  by  no  means  sufficient,  that  the  denser  material  have 
the  greater  Impedance. 
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